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Abstract 

The optimization of entanglement evolution for two-mode Gaussian pure states under one-side 
Gaussian map is studied. Even there isn't complete information about the one-side Gaussian noisy 
channel, one can still maximize the entanglement distribution by testing the channel with only two 
specific states. 
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Introduction. The study of properties about quantum entanglement has drawn much 
interest for a long time lHj]. Although initially quantum information processing(QIP) was 
studied with discrete quantum states, it was then extended to the continuous variable (CV) 

n 

quantum states[5|]. So far, many concepts and results with 2-level quantum systems have 
been extended to the continuous variable case with parallel results, such as the quantum 



re inseparability criteion[3], the degree of entanglement |8|, |9|, the entangle- 



10Nl2|. the entanglement sudden death 13], the characterization of Gaus- 



teleportation[6|, t 
ment purification 

sian mapsjl^], and so on. However, this does not mean all results with 2-level quantum 
systems can have parallel results for Gaussian states. 

Entanglement distribution is the first step towards many novel tasks in quantum commu- 
nication and Qipjlj]. In practice, there is no perfect channel for entanglement distribution. 
Naturally, how to maximize the entanglement after distribution is an important question 
in practical QIP. If we distribute the quantum entanglement by sending one part of the 
entangled state to a remote place through noisy channel, we can use the model of one-side 
noisy channel, or one-side map. 



Given the factorization law presented by Konrad et al 15[], such a maximization problem 
for entanglement distribution over one-side map does not exist for the 2x2 system because 
any one-side map will produce the same entanglement on the output states provided that 
the entanglement of the input pure states are same. The result has been experimentally 
tested [3] and also been extended recently. However, such a factorization does not hold 
for the continuous variable state as shown below. In this work, we consider the following 
problem: Initially we have a bipartite Gaussian pure state. Given a one-side Gaussian map 
(or a one-side Gaussian noisy channel), how to maximize the entanglement of the output 
state by taking a Gaussian unitary transformation on the input mode before it is sent to 
the noisy channel. We find that by testing the channel with only two different states, if a 
certain result is verified, then we can find the right Gaussian unitary transformation which 
optimizes the entanglement evolution for any input Gaussian pure state. That is to say, we 
can maximize the output entanglement even though we don't have the full information of the 
one-side map. In what follows we shall first sho w by specific example that the factorization 
law for 2 x 2 system presented by Konrad et al 15|] does not hold for Gaussian states. We 
then present an upper bound of the entanglement evolution for initial Gaussian pure states. 
Based on this, we study how to optimize the entanglement evolution over one-side Gaussian 



2 



map by taking a local Gaussian unitary transformation to the mode before sent to the noisy 
channel. 

Output entanglement of one-side Gaussian map and single-mode squeezing. Most gener- 
ally, a two-mode Gaussian pure state is 

\g(U,V,q)) = U®V\ X {q)) (1) 

and \x(q)) = \A - q 2 e qa ^ |00) (-1 < q < 1) is a two-mode squeezed state (TMSS). We 
define map $ as a Gaussian map which acts on one mode of the state only. A Gaussian map 
changes a Gaussian state to a Gaussian state only. In whatever reasonable entanglement 
measure, the entanglement of a Gaussian pure state in the form of Eq.flTJ) is uniquely deter- 
mined by q. Therefore, we define the characteristic value of entanglement of the Gaussian 
pure state p(q) = \g(U, V, q))(g(U, V, q)\ as 

E[ P {q)] = \q\ 2 - (2) 

On the other hand, any bipartite Gaussian pure state is fully characterized by its covariance 
matrix (CM). Suppose the CM of state U ® V\x(q)) is 



(3) 



\q\ 2 is uniquely determined by \A\ (the determinant of the matrix A). So, to compare 
the entanglement of two Gaussian pure state, we only need to compare \A\ value of their 
covariance matrices. 

We start with the projection operator Tk(q a ) which acts on mode k only: 

q a \n)(n\ = q t 

n=0 

This operator has an important mathematical property 




f k {q a ) = Y,1 n >){n\= q y k . (4) 



Tk(q a )(al, a k )f k 1 (q a ) = (q a a{, a k /q a ) (5) 

which shall be used latter in this paper. For simplicity, we sometimes omit the subscripts of 
states and/or operators provided that the omission does not affect the clarity. 

Define the one-mode squeezed operator S{r) = e r( - at _a2 - ) where r is a real number and 
bipartite state \i/j r (qo)) = I ® <S( r )lx(<?o))- We have 



Theorem 1. Consider the one-side map / <S> T(qi) acting on the initial state \^ r {%))- The 
entanglement for the outcome state I ® T(qi)\ijj r (qo)) is a descending function of |r|. Math- 
ematically, it is to say that if |ri| > |r 2 | then 

E[I® f ( qi )\AM)) <E[I®f( qi )\i, r2 {q ))). (6) 

This theorem actually shows that there isn't a factorization law similar to that in 2 x 2 
states for the continuous variable states, in whatever good entanglement measure. Using 
Backer- Compbell-Horsdorff (BCH) formula, up to a normalization factor, we have 

|^ r (gb)) = e -5 at i 29 o tanh(2r)+ ^ a 2 2tanh ( 2r ) + ^tA|00). (7) 

Detailed derivation of this identity is given in the appendix. Based on Eq.(jH), the one-side 
map I (g> T(qi) changes state |t/v(<7o)) into 

= e /i4 2 +/24 2 +/3«14|oo) (8) 

where fi = — |g 2 tanh(2r), f 2 = |g 2 tanh(2r), and = co ^2r) • Here we have omitted the 
normalization factor. Since we only need the covariance matrix of state \ip'), the normaliza- 
tion can be disregarded because it does not change the covariance matrix. The characteristic 
function of state p' = has the form 

C(e*i, a 2 ) = titfDiMbtM] = e~^ AdlT (9) 

where D k (ak) = e aka k~ a k ak and a = (x\, yi, x 2 , y 2 ) with a fc = -^{xk + Wk)- Writing A here 
in the form of Eq.(j3]), we find A = diag[&i, 62]) C = diag[ci,c 2 ] and B = diag[di,d 2 ] with 

U __ 1 _| 1+2/2 I. = _1_| 1-2/2 J =_Ij l + 2/l J = 

1 2 I" l+2/ 1 +2/ 2 +4/ 1 /2-/|' ° 2 2 l-2/i-2/2+4/i/2-/|' U1 2 l+2/ 1+ 2/ 2 +4/i/2-/| ' " 2 

_ 2 + 1-2/1-2/2+4/1/2-/^' Cl = I+2/1+2/2+4/1/2-/3 2 ' ° 2 = l-2/i-2/ 2 +4/i/2-/|- entan g lement 

in whatever measure of state is a rising functional of |A| and 

14 = 7 + — • (10) 

This is obviously a descending functional of |r|. 

Upper bound of entanglement evolution. Since U <8> / and J ® $ commute, the unitary 
operator £/ places no role in the entanglement evolution under one-side map I<g>$, and hence 
we only need consider the initial state \g(I,V,q)) = I <g) V\x(q)) = \<p(q)). We also define 
p G (q a ) =I®$(\<p(q a )){<p(q a )\). 



Using Eq.flSD, one easily finds \<p(q = q a qb)) = T(q a ) (g) I\(p(q b )). Since the operator 
T(q a ) <S> I and the map / (g) $ commute, there is: 

P G (q = q a .q b ) = f(q a ) <g> Ip G {q b )f\q a ) ® I. (11) 

Using entanglement of formation j^l. fl^ . we can calculate the entanglement of the state of a 



Gaussian state through its optimal decomposition form 
optimal decomposition [9 1 : 



9j. Suppose p (qb) has the following 



p G (q b ) = U 1 ® U 2 p s (q )Ul <g> U\ (12) 
Here Ui,U 2 are two local Gaussian unitaries and p s is in the form 

P s (q )= [ d 2 /3id 2 /3 2 P(/3!,/3 2 ) 



(13) 

£>(/3i,/3 2 )| X (go))(x(?o)P t (/3i,/3 2 ), 

where P(/?i, (3 2 ) is positive definite, D(^ 1 , (3 2 ) = Di(Pi) <S> D 2 ((3 2 ) is a displacement operator 
defined as Dk(flk) — e^ ka ^~^ ak . According to the definition of optimal decomposition jol, 
18], there don't exist any other t/i,[/ 2 and positive definite functional P(fii,f3 2 ) which can 
decompose p G {qb) in the form of Eq. fll2l) with a smaller \q \. The entanglement of p G {qb) 
is equal to that of a TMSS |x(<?o)) ; i- e - <7o- For the Gaussian state p G (qb) with its optimal 
decomposition of Eq.( JT2l) . we define the characteristic value of entanglement of p G (qb) as 
E[p G {q b )] = M 2 . 

Lemma 1. For any local Gaussian unitary U and operator T(q a ), we can find 8,8' and /3" 
satisfying 

T{q a )U\ <S> U 2 ■ D((3i, /3 2 )|x( ( ?o)) / x 

( 14 ) 

=U(8') <g> K{8) ■ D(/3[, (3' 2 ) ■ T(q a )S(r) <g> U 2 \ X (q Q )), 

where, S(r) is a squeezing operator defined earlier, 71(8) is a rotation operator defined by 
lZ(8)(a\a)7V {8) = {e~ %e a* , e %e a) , f3[, (3' 2 and (3i, (3 2 are related by a certain linear transfor- 
mation. 

Proof: Any local Gaussian unitary operator U\ can be decomposed into the product form 
of K{8')S{r)K{8). Also, S{r)K{8) <g> U 2 ■ D^fo) = DiP'LPZ) ■ S[r)TZ(8) <g> U 2 . Define 
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d = f(q a ) <g> / • D((3l ■ f-\q a ) <g> J, we have 



f{q a )K{6')S{r)n{6) ® I ■ /3 2 )|x(<Zo)> 

ft(0') (8) / • d ■ f(q a )S(r)n(6) ® /| x (g )) 

ft(0') (8) ft(0) • • f ( 9tt )5(r) (8) J|x(?o)). 



This completes the proof of Eg. (1141) . In the second equality above, we have used the fact 
T(q a ) and 1Z{6') commute. Also, d there is not unitary. However, using BCH formula and 
the vacuum state property ak\00) = 0, we can always construct a unitary operator D(f3[, j3' 2 ) 
so that the final equality above holds. Here /?(, /3' 2 are certain linear functions of /3 2 . 
Using Eq.(HI]) and Eq.([T2]) with Eq.flU]) we have 



£[p G (g = g a g b )] 

(8) ?7 2 • T(? a )f/i <8> Ip s U\f\q a ) ®I-I®U, 



-E 



(15) 



n{9[)®u 2 n{9 1 ) (J dPp^fcPiPufc) 

/3 2 )-f(g a ) < S(r 1 )®/|x(?o))(x(?o)|<S t (r 1 )ft(g a ) 
rf 2 /3 1 rf 2 /3 2 P(/3 1 , • T( 9a ) ® J 

lx(go))(x(go)|r t (ga)®/-£> t (/3i,/? 2 ) 

<|g a go| 2 = ^[|x(ga))(x(ga)|]-P[p G (g fe )]. 

In the third step above we have used theorem 1 for the inequality sign. This gives rise to 
the second theorem: 

Theorem 2. Using the entanglement formation as the entanglement measure, if the entan- 
glement of p G (qb) is equal to that of TMSS |x(<Zo))> the entanglement of p G (q = q a qb) must 
be not larger than that of TMSS |x(<7a<7o))- Mathematically, it is to say that if \q\ < \q b \ < 1 
we have 

E[Z®$(|y( g ))(y( g )|)] EMq))(^q)\] 
E[I®$Mq b ))&(q b )\)] ~ EMtoWfaW 
Here \f(q)) = I ® V\x(q)) as defined earlier, V can be any Gaussian unitary operator. 
Definitely, the inequality also holds if we replace \<p(q)) by \g(U,V,q)) and replace \<p(qb)) 
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by \g(U', V, qb)), and U, U' can be arbitrary unitary operators. Theorem 2 also gives rise to 
the following corollary. 

Corollary 1. Given the one-side Gaussian map I ® $, if the equality sign holds in formula 

( TTTjj) for two specific values q, qb and < \q\ < \q b \ < 1, then the equality sign there holds 

even q, q b there are replaced by any q', q", respectively, as long as \q'\, \q"\ G [\q\, 1]. 

Proof. For simplicity, we first consider the case where q is replaced by any q'. (1) suppose 

W\ e [\q\i The left side of formula ( 1T6|) is equivalent to w' -z', and w' = 

and z' = • The right side of formula ( IT6|) is equivalent to w ■ z and w = 

^n^rrm and z = ^M^WrM • Theorem 2 itself says that w' < w and z' < z. If the 

E[\<p(q')){<pW)\] E[\<P(<lb)){<p(qb)\] 

equality sign holds in formula (116H . we have w' ■ z' — w ■ z hence we must have w = w' and 
z = z' which is just corollary 1 in the case q is replaced by q'. (2) Suppose \q'\ > \qb\. As we 
have already known, p G (q) = T(q a ) (g) Ip G (qb). Consider Eq. (|14p . Unitary U\ in the optimal 
decomposition of Eq. (|T2l must be a rotation operator only, i.e., it contains no squeezing, 
for, otherwise, according to theorem 1, E(p > ((/)) i s strictly less than q\q\ which means the 
equality in formula ( |T6l) does not hold. 

We denote q' = qb/q c and \q c \ < 1. We have 

p G (q = qb/qc) 
= f- l {q c )®I P G {q b ) ff- 1 ^)® J Xt 



= f~ x {q c ) U 2P S U\ ® U\ ■ f- x {q c ) ® I 

= H 1 ®U 2 - j d 2 (3 1 d 2 (3 2 P{(3 1 ,!3 2 )b{(3[,(3' 2 ) 

\x(qo/qc))(x(qo/qc)\DH(3[, &) • ^ ® t4- (17) 

Here we have used T _1 (g c ) g) = 9'9c)) = |x(9'))- ^ e have used the optimal decompo- 

sition for p G (qb) in the second equality, and lemma 1 in the last equality above. Eq. fllTl) is 
one possible decomposition of the state p G {q'), but not necessarily the optimized decompo- 
sition. Therefore, E[p G (q' = qb/qc)} < |<Zo| 2 /|9c| 2 = |<7'| 2 /|9&| 2 " E[p (qb)]- On the other hand, 
according to theorem 2, we further obtain that E[p G (q b = q'q c )} < |%| 2 /k'| 2 " E[p G (q')j. 
Remark: Since here \q'\ > g&, sign < should be replaced by sign > in formula f|T6|) . when q 
is replaced by q'. These two inequalities and result of (1) lead to 

E[p G (q')} E[\x(q'))(x(q')\] 



E[p G (q b )] E[\ X (q b ))( X (q b )\Y 



(18) 
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for any q' provided that \q\ < \q'\ < 1. Replacing symbol q' above by symbol q", we have 
another equation. Comparing these two equations we conclude corollary 1. 
Lemma 2: Given any Gaussian unitaries U, V, we have 

E[I®${U®V\<f> + ){<j> + \tf ®V^)] = £[/<g>$(|0 + ))]. (19) 

Here \<fi + ) is the maximally entangled state defined as the simultaneous eigenstate of position 
difference X\ — x<i and momentum sum pi + p 2 , with both eigenvalues being 0. Also, when 
q = 1, the state \x(q)) — We shall use the following fact. 

Fact 1: For any local Gaussian unitary operators U and V, we can always find another 
Gaussian unitary operator V so that 

U ®V\<P + ) = V ®I\<p + ). (20) 

Proof: Any local Gaussian unitary operator can be decomposed into the product form of 

TZ{e')S{r)n{e). For any TMSS \x(q)} we have H(6i) ® H(9 2 )\x(q)) = I ® K(0i + 

For the maximally TMSS \<ft + ) we have S(r) (8> S(r)\<f) + ) = \<fr + ), for, the both sides are the 

simultaneous eigenstates of position difference and momentum sum, with both eigenvalues 

being 0. This also means S(r) <8> = I (g) 5 t (r)|0 + ). Suppose V = n(6' B )S{r B )n(6 B ), 

then 

U ®V\<P + ) = V®I\<p + ) (21) 

where V = U1Z{6 B )S\r B )1Z{6' B ) . This completes the proof of Eq. (l20|) . If the equality sign 
in formula ffTBT) holds, we can apply corollary 1 of theorem 2 through replacing q^ by 1 and 
we obtain that E[p (q')\ = \q'\ 2 ■ E[I ® $(|0 + ))]. On the other hand, by using theorem 2 and 
lemma 2 we have E[p G (q')] < \q'\ 2 ■ E[I <S> $(|0 + ))]. This means 

E[p G (q')} = max{E[J ® $(|^(J, V, q')))]} (22) 

where p G (q') = I <S>$(\g(I ,V, q 1 )) (g(I ,V, q')\) as defined earlier, {V^'} is the set containing all 
single-mode Gaussian unitary transformations. The equality holds for any q' provided that 
the equality of formula ffTBT) holds for two specific values q, q^ and \q'\ > \q\. We arrive at 
the following major conclusion of this Letter: 

Major conclusion: Suppose that we have a TMSS \xW))- We want to maximize the entan- 
glement distribution over a one-side Gaussian map I <E> $ by taking local Gaussian unitary 

8 



operation I £g> V before entanglement distribution. Although we don't have complete in- 
formation of the map / ® $, it's still possible for us to find out a specific Gaussian unitary 
operation V so that the entanglement distribution is maximized over all V, for an initial 
state \xW)) with any \q'\ > \q\, as long as we can find two specific values \q^\ > \q\, such 
that the equality sign in formula ( II 6p holds. Obviously, the conclusion is also correct for 
any initial state which is a Gaussian pure state. 

The conclusion actually says that, in verifying that V can maximize the entanglement 
distribution for all initial states > |?|}> we om y need to verify the equality sign 

of formula (1161) for two specific values. 

Experimental proposal. To experimentally test our major conclusion, we can consider the 
following beamsplitter channel: Initially, beams 1 and 2 are in a TMSS, which is the initial 
bipartite Gaussian pure state. Beam 3 is in a squeezed thermal state p3 = S(u3)p t hS^ (u 3 ) 
here S(u) is a squeezing operator defined by S(u)(x ,p)S^(u) = (ux,p/u) and p th is a ther- 
mal state whose CM is diag[&3,&3]. Beam 3 together with the beamsplitter makes the 
one-side Gaussian channel. A beamsplitter will transform £2,^3 by Ub(%2, x^jU^ 1 — > 

(cos 9 sin 9 \ 
. In an experiment, we can take, e.g., q = 0.02 and q^ = 0.5, test- 
— sin 9 cos 9 I 

ing with many different V we should find that the equality sign in formula ( JTTjj) can hold 
with V = S(v,2 = M3) • Our major conclusion is verified if we can find that the same 
V = S(us) always maximizes the output entanglement for any input state |x(</)) provided 
that \q'\ > 0.02. Numerical calculation is shown in the following figure. 

The numerical example with q=2/3, G=jt/6, b3=1, u3=3. 

0.27 

0.26 

_ 0.25 
s= 

E 

1 0.24 
01 

1 0.23 
H 0.22 
0.21 
0.2 

1 23456789 10 
The squeezing factor u2. 

FIG. 1: The entanglement with different squeezing factor ui- The maximum entanglement obtained 
when U2 = u% = 3. Here we set U3 = 3 and q' = 2/3,8 = tt/6, 63 = 1. 
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In summary, we present an upper bound of the entanglement evolution of a 2-mode 
Gaussian pure state under one-side Gaussian map. We show that one can maximize the 
entanglement distribution over an unknown one-side Gaussian noisy channel by testing the 
channel with only two specific states. An experimental scheme is proposed. 
Acknowledgement. This work was supported in part by the National Basic Research Pro- 
gram of China grant nos 2007CB907900 and 2007CB807901, NSFC grant number 60725416, 
and China Hi-Tech program grant no. 2006AA01Z420. 

Appendix. Details of the proof of Eq.^jJ. We will use the following lemma. 
Lemma 2. If A and B are two noncommuting operators that satisfy the conditions 



[A,[A,B}} = [B, [A,B}] = 0, 



(23) 



then 



D A+B = e A e B e -\[A,B] 



This is a special case of the Baker-Hausdorff theorem of group theory 1^1 . 
The squeezing operator S(r) = e r( - at ~ a2s> can be normally ordered as 20] 

,t 2 



S(r) 



1 



a/ cosh(2r) 



cxp 



tanh(2r) 



■-a 2 tanh(2r) 



• exp [— a^a(ln(cosh(2r)))] exp 
We neglect the constant of normalization in all the following calculation. 

I ® S(r)\x(qo)) - 



K4 -a2) e <2oal<4|(]0) 



loo) 



^q a\(a\ cosh(2r) — a 2 smh{2r)) ^r(a^ — a|) | 

^q a\(a\ cosh(2r)-a 2 sinh(2r)) e ia^ 2 tanh(2r) | qq^ 
I„t 2 

e 2 



a| tanh(2r) ggo a I { a 2 cosh(2r) — [a2+a\ tanh(2r)] sinh(2r)} |QQ^ 
\af tanh(2r) e go4( C o 3 h(2r) - a2 sinh ( 2 ^)) | qq^ 



|4" tanh(2r) £ Z^\t.t) e ~^\'<^ tanh(2r) | gg^ 



1o44 l„t 2 2 ( 



(24) 



(25) 



(26) 



This is just Eq.(H). In the last equality we have used lemma 2. This completes the proof of 
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